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Abstract 

The exactly solvable scalar-tensor potential of the four-component Dirac 
equation has been obtained by the Darboux transformation method. The 
constructed potential has been interpreted in terms of nucleon-nuclear and 
Schwinger interactions of neutral particles with lattice sites during their chan- 
neling in the nonmagnetic crystal. The family of exactly solvable interaction 
Hamiltonians of a Schwinger type is obtained by means of the Darboux trans- 
formation chain. The analytic structure of the Lyapunov function of periodic 
continuation for each of the Hamiltonians of the family is considered. 



1. Introduction 

The Darboux (Moutard) transformation [H [2j |3] is widely used for constructing solv- 
able quantum mechanic models [H 0, El [7j . It is known as the method of construction 
of the solvable potentials for the one- dimensional Schrodinger equation P, [9j QUI E] 
and one for two-component case was considered in [12]. The quaternionic factoriza- 
tion of the Schrodinger equation was obtained in [13]. 

As was shown in [TH [15l UHl [T71 HHJ H5] . the exactly solvable models for the 
one-dimensional Dirac equation with the Lorentz scalar potential employ the close 
relation between this equation and a pair of the Schrodinger equations. The solution 
of these equations, in turn, is based on well established algebraic techniques. In 
papers J2H1 12H E21 [23] the Darboux transformation was applied directly to the one- 
dimensional two-component Dirac equation, where as for applications it is often 
important the four-component Dirac equation [2ll [25] . 

In works [261 EZ] the attempts to use the Darboux transformation to the four- 
component Dirac equation have been done. However, as a result, the Darboux 
transformation was applied to two-component equations. 

In [13] the Darboux transformation of the Dirac equation was reduced to the 
quaternionic factorization of the Schrodinger equation. The Darboux — Crum method 
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[H [28] was applied to the Schrddinger equation in [29] and to the two-component 
Dirac equation in [2TJ . 

In the present work we apply the Darboux transformation to the one-dimensional 
four-component Dirac equation, thus we make contribution for applying the Dar- 
boux transformation to the full Dirac equation. Also we construct a chain of the 
Darboux — Crum transformations for the four- component Dirac equation. We try 
to generalize these results to construct an exactly solvable interaction Hamiltonian, 
which allows one to solve analytically the one- dimensional four- component Dirac 
equation. 

The structure of the present paper is as follows. In section 2 we construct a Dar- 
boux transformation for the free Hamiltonian of the one- dimensional four- component 
Dirac equation and generate an exactly solvable transparent scalar-tensor potential. 
We show that the tensor part of the external potential has the meaning of the elec- 
trical strength. In section 3 we consider a some special choice of transformation 
parameters. In section 4 we construct a chain of the Darboux — Crum transforma- 
tions for the four-component Dirac equation. In section 5 we briefly discuss the 
structure of the Lyapunov functions for the problem of the interaction with the 
periodic potentials. Finally, in section 6 we interpret the interaction of a spin-1/2 
neutral massive particle with an external electrical field through its magnetic mo- 
ment as coherent Schwinger interaction and try to apply these results to neutral 
particles channeling in crystals. 



2. Darboux transformation of the one-dimensional 
four-component Dirac equation 

The stationary one-dimensional Dirac equation reads: 

Htp — Eip, H = H + Hj, (1) 
H = —ia\d x + /?m, Hj = H[(x), 

where ai = ^ ® °^ j , (3 = ^ ^ 2 ^ j , 1 2 = ^ J ^ j ; <J\ is the Pauli matrix; 

Ho is the free Dirac Hamiltonian [261 E3 [3Ql EI] and Hi is the part of the Hamiltonian 
describing an interaction of the spin-1/2 particle with some external field. The 
solution of the equation ([1]) in the absence of the interaction is simple and well 
known for any E. 

Let us consider the Darboux transformation 

HL = LH , ?p = Lip , Hip = Eip (2) 

with the operator L of the form 

L = d x — UxUT 1 . (3) 

Here 4 x 4-matrix u obeys the equation 

H Q u = uA, (4) 
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the A is any nonsingular matrix, u x = du/dx. 

Such transformation generates the interaction Hamiltonian 

Hj = H — H = —i(ctiu x u~ l — u x u~ a\). (5) 

If we verify that the function is the solution of Eq. where Hj is defined by 

(EJ), we see hat u is definite by (TJ]). 

We choose A to be diagonal matrix of the form 

A=[A 1 (I + /?) + A 2 (I-/3)]/2, Ai, 2 <m. (6) 

1 2 



Here / - - , 

U I2 

Let us construct the matrix u: 

u = al + bai + c(3 + dj, 7 = (7) 

where 

a = (/ii + /i 3 )/2, 6 = (/i 2 + /i 4 )/2, c = (/ii - /i 3 )/2, d = (// 4 - // 2 )/2, (8) 

/ii = cosh (fcix + 0i), yU 3 = cosh (fc 2 x + 2 ), (9) 

-ik 2 sinh (fc 2 x + 2 ) -i/ci sinh (fcix + 0i) 

^2 = : , /i4 = , (10) 

A 2 — m Ai + m 



ki : 2 — ±y?ri 2 — A^ 2 ; 0i )2 are arbitrary complex numbers. It is easy to check that 

these matrix is solution of the Eq. (00). 
The inverse matrix reads: 

iT 1 = [al - bat -cP- dj]/D, D = a 2 — b 2 — c 2 + d 2 . (11) 

With the help of simple calculations it can be proved that 

Hj = H s + H T , H s = (3Vs, H T = i 1 V T . (12) 

Here Vg, Vr are functions of x. The meaning of Vs is transparent. It is an external 
scalar potential. The tensor part of the interaction Hamiltonian can be consid- 
ered as particular case (A^ = (A (xi),0)) of the more general tensor Hamiltonian 

H T = gho^^F^, (13) 
<V = lulv-lul^ (14) 
F, v = (15) 

where 7 M are the Dirac matrices, A^ is an external electromagnetic field and g is a 
coupling constant proportional to magnetic moment of the particle. 

Such Hamiltonian describes the interaction of a spin- 1/2 neutral massive particle 
(e.g. neutron) with an external electrostatic field. 

From above consideration we conclude that the tensor potential Vr has meaning 
of the electrical strength 

(9/4 

V T ~ E x = (16) 
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Figure 1: The scalar part of the external potential Vj(x) = Vs(x) + Vr(x): Vs(x) 
-2k 2 /(m + \cosh(2kx + 2<p)), k = ^m 2 - A 2 , m = 1, A = 1/2. 

3. Some special choice of parameters 

Case 1. Let Ai = A 2 = A, 4>i — 02 = 0- In this case Vt = 

2A; 2 



m + A cosh (2kx + 2(f)) 



k = \J m 2 — A 2 . 



(17) 



This result coincides with the result of [32] at = [In (m + k)/(m — k)]/4. 
The function ^(x) is shown for fixed values m — 1 and A = 1/2 in Figure [U 
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Figure 2: The tensor part of the Vi(x): Vr(x) = ki tanh(A; 2 x), k 2 = \Jm 2 — A 
m = 1, A 2 = 1/2. 



Case 2. Let Ai = m, 0i = 2 = 0. In this case 

Vs = — A 2 — m = const 
which leads to the redefinition of the particle mass 

m — > — A 2 . 



;is) 



(19) 
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The tensor potential (electrical strength) is 

Vt(x) = &2 tanh (k%x) ~ E x . 



(20) 



A pictorial representation of the potential Vt{x) for fixed values m = 1 and A2 = 1/2 
is given in Figure [2j 

Such electrical field is produced by the charge distribution 

" w ~ (21) 

The charge distribution p for values m = 1 and A2 = 1/2 is shown in Figure [3J 
It is easy to check that both Vs and Vr are transparent potentials. 

4. Simple chain of the Darboux — Crum transfor- 
mations 



As the next step consider a chain of the Darboux — Crum transformations for the 
four-component Dirac equation ([T]). Let us put in expressions ©, ffTUl) 4>i = 02 = 0, 
fci = 0, k 2 = k, A2 = £1 = \fmF^W . This leads to the Hamiltonian of the 
following form H\ = —ia\d x — S\f3 + ry tanh (fcr). 
Consider the matrix 



where 



u 



(i) = a «I + b^at + c«/3 + dW T , 



- ^)/2, 



cosh (fcx), 
cosh 2 (for) 



f. ( " = t/4" + f i")A 

d« = (rf' - ^")/2, 
— 3i/c sinh (fcx) cosh (/ex) 



£2 - £\ 



0. 



e 2 = Vm 2 — Ak 2 . 



(22) 
(23) 

(24) 
(25) 



p(x) 
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Figure 3: The charge distribution p(x) = 1/ cosh 2 (fc 2 ^), &2 = \J m 2 — A|, m 
A 2 = 1/2. 
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The matrix (12211 obey the equation 



HiUi = u\Ai, Ai 



1 



(/ + /3)aW + (/-/3)aW 



r A (1) - f 



A 1 — — ci, /v 2 

Consider the intertwining relation 

Ll#l = i?2^l, 

where 

L = d x -u<puW-\ 

Then 

H 2 — Hi = —i[ai,u^uP~) 1 ] = (3(e\ — £2) + 27&tanh (kx) 
from which it follows 

H 2 = —iot\d x — s^ft + ry2A;tanh (foe). 



(26) 
(27) 

(28) 

(29) 

(30) 

(31) 




Figure 4: The tensor part of the interaction Hamiltonian ( 1321) : Vr(x) = nk tanh(fcx) 
k = \J m 2 — Xf, m = 1, A2 = 1/2 (line 1: n = 1; line 2: n = 3; line 3: n = 5) 



Let 



= —ia.\d x — e n j3 + rynfctanh (kx), e n = \lm 2 — n 2 k 2 



and 



where 



a (n) = (//i n) +^ n) )/2, 6 (B) = (A4 B) +A*i B) )/2, 
= (/4 n) - /4 n) )A rfW = ( / if- /U r ) )/2, 



(32) 
(33) 
(34) 
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(n) ,n/, \ (i) -i{2n + l)fcsmh(fcr)cosh (n 1} (kx) , s 

fi\ = cosh {kx), fi 2 = , (35) 

^ = cosh (n+1) {kx), = 0. (36) 

Then u n obey the equation 

\(I + + (I - P)>£>] , (37) 



H n u n u n \ n , A n ^ 



A! n) = s n , \ {n) =e n+l . (38) 
The solutions of the equation 

H n+1 i> s = E s ij s (39) 
corresponding the discrete levels are of the following form: 



si 



{<Psl,<frs2,(fis3,<fis4,y, 

AP^{tanh (kx)}, s2 = PP^{tanh (kx)}, 
fi s3 = OP* +1 {tanh(fcx)}, (j) s4 = PP* +1 {tanh {kx)}, 
A B _ k{n + s + 1) 
D C e n+ i + E s 



(40) 



E s = sign{s)e s , e s = \Jm 2 — s 2 k 2 , (41) 
s = -n, -Ti + l, ...,-2,-1,1,2, ...,n,n + l, (42) 

where P^ are associated Legendre polynomials [33], P and C are arbitrary nonzero 
numbers. 

The solutions of the equation 

H n+1 ip = Eip, \E\ > m, (43) 

corresponding to the continuous spectrum can be obtained from (HOI) by the substi- 
tution 



P* (n+1) {tanh(fcr)} =^ P£ (n+1) {tanh (fcz)}, /i = ±ip/k, p = VE 2 - m 2 , (44) 



P<^> {ta h (fa)} = • F (-»; n + 1; 1 T f ; i^M) . (45) 

From the last expression it is evident that "potentials" of the form — /3e n +i r yn tanh {kx) 
are transparent. 



5. Exactly solvable periodic potentials 

Starting with exactly solvable potential for full real axis one can construct the exactly 
solvable periodic potentials in the manner described in the papers [TQl 1221 EI] ■ 
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It is well known that the continuous spectrum of such potentials have a band 
structure. For the investigation of this structure it is necessary to construct so-called 
Lyapunov function D(E) (see for details pHl E2, EI])- 

The values of energy E which obey the equation \D(E)\ > 2 belong to the 
forbidden bands. All others belong to the allowed bands. 

Omitting simple but cumbersome calculations we present here only final result 
for the some lowest values of the n: 



D n (E) = 2[cos(2pa)A n (p, tx) - (ktx/p) sm(2pa)}B n (p, tx), (46) 

t x = taoh(A;a), a = T/2, (47) 

where T is the parameter of periodization [TDl [221 Elj- In general case (of the 
arbitrary n) 

A ^ = Im^ n {uM) ^ = M„.(Mi) 

V t\ 

$ n = (n - iv) ■ F (-n + 1, n; 1 - iv\ ^-^J ■ F ^-n, n + l;l + iv; , (49) 

v — p/k. 

In the case ka >> 1 (1 — tx << 1) 1), it is readily found from (??) and (??) 




Figure 5: Lyapunov functions D nj n — 1, 3, with ka = 6. 



D = 2 cos(0„. -arctan(nfc/p)) 
cos(arctan(n/c/p)) ' 

n 

n = 2[pq + arctan(sfc/p)]. (51) 

The solution of inequality \D n (E)\ > 2 is trivial in this case. 

Lyapunov functions D n (E) with n — 1, 3 are presented in Fig. [5] for fca = 6. 
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6. Discussion 



We show that the Darboux transformation for the free one-dimensional four-compo- 
nent Dirac Hamiltonian generates the interaction that can takes place in the real 
world. 

The scalar potential arises in the theory of nuclear-nuclear interactions [351 136] 
when sigma-exchange is taken into account [371 138] and can be considered as a part 
of the nuclear potential in the channeling theory [391 HO] ■ 

The interaction of a neutral massive particle with the spin 1/2 (e.g. neutron) 
with an external electrostatic field generated by the Darboux transformation can 
be interpreted as Schwinger interaction of neutron with atoms of the lattice in its 
channeling in the nonmagnetic crystals [391 EJ H21 H31 0U 051 HE] . 

The origin of the Schwinger interaction is due to relativistic effect of appear- 
ance of the lattice magnetic field in the neutron rest frame, when it passes through 
the electrical field of the crystal. Due to proper magnetic moment of neutron the 
appearance of such field leads to magnetic interaction, which allows to order and 
orientate movement of neutron as in the case of channeling in nonmagnetic crystal 

[391 SQl SZl SHI SS]- 

Really, in the coordinate system connected with the neutron moving in the elec- 
trical field of the crystal the magnetic field H = [E,v\ appears. In the case of the 
plane channel the value of the E averaged by the area of elementary cell goes to 
zero. As a result, the mean value of the H obtains its maximum near the plane. In 
the middle of two planes H goes to zero, changes the sign and increases to maximum 
value near the next plane of the channel (see Figures 2 and 4). 

Moving across the plane the effective field H, due to the changing of the sign of 
the electrical field vector E, change the sign and its space behaviour repeating the 
behaviour in the previous channel (Figure 5). 

V T (x) 




Figure 6: Periodic continuation of the potential ( J20l) for k 2 = 1. 

In accordance with the accepted terminology, the considered neutron-cell interac- 
tion may be characterized as coherent Schwinger interaction and is the generalization 
of well known neutron-atom Schwinger interaction |39j . 
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